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Abstract: The fc-string tensions are explored in the 4-d [U(l)]^ _1 -invariant dual Abelian- 
Higgs-type theory. In the London limit of this theory, the Casimir scaling is found in the 
approximation when small-sized closed dual strings are disregarded. When these strings are 
treated in the dilute-plasma approximation, explicit corrections to the Casimir scaling are 
found. The leading correction due to the deviation from the London limit is also derived. 
Its iV-ality dependence turns out to be the same as that of the first non-trivial correction 
produced by closed strings. It also turns out that this iV-ality dependence coincides with 
that of the leading correction to the fc-string tension, which emerges by way of the non- 
diluteness of the monopole plasma in the 3-d SU(iV) Georgi-Glashow model. Finally, we 
prove that, in the latter model, Casimir scaling holds even at monopole densities close to 
the mean one, provided the string world sheet is flat. 
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1. Introduction 

The spectrum of fc-string tensions has been extensively explored in recent years, both on 
the lattice ]l]]-[j7|] and in the continuum limit [p[-Jl3|]. A fc-string is a string joining k quarks 
with k antiquarks. Alternatively, it can be defined as a string between sources that carry 
charge k (iV-ality) with respect to the center group, Zjy, of the original SU(iV) group. The 
knowledge of the tensions of A;-strings, a^, is likely to shed some light on the dynamics of 
confinement by providing some new "phenomenological" input that models of confinement 
have to satisfy, see e.g. [l5[| . 

Some properties of the /c-string spectrum follow directly from first principles. The 
ratio of to the tension of the fundamental string, a%, in an SU(iV) pure gauge theory is 
invariant under the interchange of quarks and antiquarks, k *-* (N — k). Moreover, in the 
large- iV limit, the interactions between fundamental strings inside the composite /c-string 
are suppressed, it is therefore expected that in this limit their spectrum fulfils the condition 



k. Finally, the corrections to this large-A limit are expected to appear in powers 



of 1/N 2 I 



Analytical calculations in supersymmetric theories [16, 11], or in some approximations 



to Yang-Mills theories based on the string/gauge duality [17, |jq| , yield the so-called sine 
scaling for the fc-string spectrum: — = s ™( fc7r /^0 i n t wo dimensions, one can show that the 

o\ sm(7r/iV) ' 

string tension ratio obeys an exact Casimir scaling: ^ = k ^_^ ■ The same result is ob- 
tained for the four-dimensional pure gauge theory, when the computations are based either 
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on dimensional reduction [19| or on the low orders of some perturbative expansion; see e.g. 
Ref. H for the fc-string spectrum in the strong coupling Hamiltonian formulation. Casimir 
scaling also appears in some models of the QCD vacuum, such as the stochastic vacuum 
model of QCD ||, || and the [U^l)]^" 1 gauge-invariant Abelian-projected theory [[!(]]. In 
3-d, Casimir scaling with certain corrections has been found in Ref. [[l3|]. 

It is important to stress that no analytical computation of the fc-string spectrum has 
been done so far that could be directly applied to the case of pure gauge theory at weak 
coupling; neither the sine nor the Casimir scaling can be considered as exact results in non- 
supersymmetric Yang-Mills theories. While both formulae yield the correct limit as N — > 
oo at fixed k; at subleading order, the sine scaling shows the expected 1/N 2 corrections, 
whereas the Casimir scaling leads to corrections in powers of 1/N, thus contradicting 
the results of Ref. [12|. Lattice results show that both formulae give an approximate 
but satisfactory description of the numerical data. Different lattice simulations display 
some discrepancy in the spectrum, which are likely to arise from systematic errors in the 
computation of the string tension. At the current level of accuracy, lattice data do not allow 
a clear-cut distinction between the two behaviors; of course there is no theoretical reason 
to expect any of the two to be exact, and it would be more interesting if lattice results 
could become accurate enough for one to be able to pinpoint the parametric behavior of 
the subleading corrections to the large- N limit. 

This paper parallels in its spirit Ref. [13], since it also deals with a confining Abelian- 
type theory, although a 4-d one. More specifically, we explore in this work /c-string tensions 
in the 4-d [L^l)]^ -1 gauge-invariant Abelian-projected theory, which is formulated in terms 
of dual magnetic Abelian gauge fields, neglecting the off-diagonal degrees of freedom. Un- 
like Ref. where this was done in the Bogomol'nyi limit and on the basis of the analysis 
of the classical string solutions, here we will consider the model in the London limit (and 
in its vicinity), which corresponds to an extreme type-II dual superconductor; the ^-string 
tensions will be derived from the string representation of the partition function. Such a 
representation means a reformulation of the theory in terms of the path integral over closed 
dual strings, which are always present in the theory and interact with the external /c-string. 
As a brief historical remark, let us mention that, for the usual (dual) Abelian Higgs model, 
the reformulation of the partition function in terms of closed strings has been performed 
by many authors (e.g. in Ref. 2C] in 4-d and in Ref. [21] in 3-d), in particular with the 
applications to the stochastic vacuum model [^, 23 1; the Jacobian of the transformation 
from field to string variables has been discussed in Ref. 24]; the SU(iV) generalization has 
been explored, in particular for studies of the 0-term [^o|] [sec Ref. for investigations of 
the 0-term in the SU(2)- and SU(3)-inspired cases] and for the purposes of further appli- 
cations to the stochastic vacuum model at N = 3 |27|, [H| [see also [28] for related studies 
at N = 3]; the corrections emerging in the vicinity of the London limit have also been 
explored |3C 



The important fact is that closed strings are short-lived (virtual) objects [31], whose 
typical sizes are much smaller than both distances between them and size of the /c-string 
world sheet. Therefore, since world-sheet tensors (also called vorticity tensor currents) of 
closed strings enter the final action in the linear combinations with the world-sheet tensor 
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of the fc-string, in the leading semi-classical approximation the interaction of closed strings 
with the fc-string can be merely disregarded. This is precisely the approximation in which 



the SU(2)- and SU(3)-inspired models have been considered in Refs. [22, 23] and 
respectively. This approximation can be further improved by treating closed strings in 
the dilute-plasma approximation |29]. In this paper, specifically in Section 5, we will also 



account for effects produced by the dilute plasma of closed dual strings. 

The sketch of the paper is as follows. In the next section, the model under study 
will be described. In Section 3, will be derived in the London limit. In Section 4, we 
will explore corrections to at, which emerge in the vicinity of the London limit, i.e. when 
Higgs bosons are not infinitely heavy. In Section 5, we will return to the London limit 
and consider another type of corrections, namely those produced by closed dual strings. In 
Section 6, the main results of the paper will be discussed once again. In Appendix A, some 
estimates related to the main part of the text will be performed. Finally, Appendix B is 
devoted to the 3-d SU(A r ) Georgi-Glashow model; it is complimentary to Ref. |l3|], where 
Casimir scaling in this model has been proved for an arbitrarily shaped surface, but at 
monopole densities much lower than the mean one. Here we prove that Casimir scaling 
holds already at the mean density, provided the surface is flat. 

2. The model 

The model we are going to deal with is the generalization of the SU(3)-inspired dual 



Abelian-Higgs-type theory [32] to the case of arbitrary N |T(], g5|. The monopole conden- 
sation is modelled in it by the assumption that monopoles form condensates of the dual 
Higgs fields. This model can naturally be called an effective [U^l)]^ -1 gauge-invariant 



Abelian-projected theory. Its Euclidean partition function reads: 
\v\$i\ve}j VB^5 \12 d ij exp |"/ d ^ x 

+ Yl - ^mq*B M ) + A (|$i| 2 - 7? 2 ) 2 




(2.1) 



Here, the index % runs from 1 to the number of positive roots of the SU(A r )-group, that 
is N(N — l)/2. Next, g m is the magnetic coupling constant related to the electric one, 
g, by means of the Dirac quantization condition g m g = Aim. In what follows, we will for 
simplicity restrict ourselves to the monopoles possessing the minimal charge only, i.e. set 
n = 1, although the generalization to an arbitrary n is straightforward. Note that the 
origin of root vectors in eq. ( |2.1[ ) is the fact that monopole charges are distributed along 
them. Further, <£j = |$j| e l6i are the dual Higgs fields, which describe the condensates of 
monopoles, and = c^Bj, — 9^B^ is the field-strength tensor of the (N — l)-component 
"magnetic" potential B^. The latter is dual to the "electric" potential, whose components 
are diagonal gluons. Since the SU(iV)-group is special, the phases 9j of the dual Higgs fields 
are related to each other by the constraint ^ 6i = 0, which is imposed by introducing the 
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corresponding <5-function into the r.h.s. of eq. (|2.1[) , Next, = \e pv \ p O\ p , and F pv is 
the field-strength tensor of k test quarks with colors dj, which move along a contour C. This 

k 

tensor obeys the equation S^F* = gMjjVj where = m ai , j p (x) = § dx^(r)5(x — 

i=l c 

x(t)), m ai is a weight vector of the group SU(iV), and Oj may take values 1, . . . , N. Thus 
F^j, = — gMfcS^, where = J do" Atl y(3;(^))(5(x — x(£)) is the vorticity tensor current 

associated with the world-sheet £ of the open electric string, bounded by the contour 
C. From now on, we will omit the normalization constant in front of all the functional 
integrals, implying for any k the normalization condition Z\. [C = 0] = 1. 

Next, the phases of the dual Higgs fields can be decomposed into multivalued and 
single- valued (also called singular and regular, respectively) parts, Oi = 8f ne + #[ eg . The 
fields 0^ mg describing closed dual strings are related to the world-sheets Sj of these strings 
by means of the equation 



e^x p d x d p ef ns (x) = 27rE* B/ (x) = 2tt J da pv 5 (x - x®(0 ) . (2.2 



This equation is the covariant formulation of the 4-d analogue of the Stokes' theorem for 
du&i, written in the local form. In eq. ( |2.2D , = x p \^) is a vector, that parameter- 

izes the world-sheet Sj with £ = (S, 1 ,^ 2 ) standing for the 2-d coordinate. As far as the 
regular parts of the phases, #[ eg , are concerned, these describe single-valued fluctuations 
around closed strings, which are described by # ? - mg . Note that, owing to the one-to-one 
correspondence between 6f ng and Sj, established by eq. ( |2.2[ ), the integration over #" ng 
is implied in the sense of a certain prescription of the summation over world-sheets of 
closed strings. One of such prescriptions, corresponding to the dilute plasma of closed 
strings, will be considered in Section 5. Further, by virtue of eq. (2.2), it is also possi- 



ble to demonstrate that the integration measure T>9i becomes factorized into the product 
V6f ng V9^ eg . Apparently, the constraint imposed by the 5-function S I ^29% ) becomes also 



sphtas5(EC S )^(EC S ) 



Let us further expand in eq. (2.1) as |$,| = r\ + ■% and perform the gauge trans 



formation <7 m qjB" ew = g m qjB /t — d p 9i, noticing that, according to eq. ( |2.2| ), (d^dy 
d u d p )9^ mg = 27rS* . The constraint ^ 9" ng = leads to the constraint for world-sheets of 



closed strings, which can be imposed by the 5-function 5 I ^2 £5^,1 . Instead, the constraint 

imposed by 5 #J cg ^ can be lifted to the exponent upon the introduction of some La- 
grange multiplier. By virtue of the fact that ^ q, = 0, the integration over this multiplier 



can be shown [27] to eventually result in an inessential overall constant factor. If we further 



use the orthonormality of roots, 
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E«?9f = T^' (2 - 3) 



2 

where a, /3 = 1, . . . , N — 1, we obtain the following Lagrangian 

2 m 2 



+ ^B 2 + 



(2.4) 



The masses of the dual vector boson and the dual Higgs field here read m = gmVv N an d 
ruH = 2r/vA- 

3. London limit 

In terms of the Landau-Ginzburg parameter k = -p^, the London limit (LL) is defined 
as ln/t 3> 1. Let us now proceed with the string representation of the model ( |2.4| ) in this 
limit. Notice that, since we would like our model to be consistent with QCD, we must have 
g = \J X/N, where A remains finite in the large- iV limit. The definition of the London limit 
then yields 

(2vreA0 2 , . 

A>^ — =^-. (3.1) 
A 

In general, in order to have confinement (i.e. the type-II dual superconductivity) in the 
large- iV limit, one should demand that A grows with N at least as 0(N 2 ). If A grows 
with N faster than iV 2 , A = 0(N 2+t ) where e > 0, k grows with N too, n = 0(iW 2 ), 
making the London limit deeper. In what follows, we will adopt the minimal requirement, 
A = 0{N 2 ), necessary for confinement in the large- N limit of the model (|2.1|). As we have 
seen, only in this case, k is N- independent. 

In the London limit, the partition function of our model has the form 



o-LL 



J (IT [<®U ) 6 fe K^j J exp (- J d A xC\ h ^ , (3.2) 
where the Lagrangian reads 



11 



+ ^B 2 . 



The symbol [cZE^,] in eq. (3.2) is a formal expression for the sum over string world-sheets, 
whose concrete form will be specified below. 

A natural way to satisfy the constraint ^ S* = is to set = qjS M ^, where S^ u 

i 

are then no longer subject to any constraint. The Lagrangian then takes the form 
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4a 2 
^ar? 



5mF^ - 4vrM fc S^ + 2^S Ml/ + — B 



m 



To perform the integration over B^, let us linearize the squares in this Lagrangian upon 
the introduction of two auxiliary fields, hau and k„, as follows: 



*-k 



1 



The integration over B„ then yields 



and we obtain 



■m 



V2 



8 h 



(3.3) 



(3.4) 



(3.5) 



To perform the integration over the so-called Kalb-Ramond field h^ u , notice that the 
general solution to eq. (|3.4|) with respect to h.^ v {x) reads 



x) 



-VZme^xp&x J d A x'D (x - x')k p (x') + d p C u (x) - d v C p (x), 



where is an arbitrary vector field, and Dq(x) = 1/(4tt x ) is the Coulomb propagator. 



To fix Cm = is equivalent to impose the constraint dnh 



- 0. On the other hand, 
because of the coupling of to the open-string world-sheet, — yM^h^S^, the saddle- 
point value of h^, which saturates the respective Gaussian integral, does not obey the 
constraint d^h^ = 0. To make the integration consistent, we should therefore promote 
h^j, by making C M non- vanishing. For such h^, (d^h^) 2 in eq. (|3.5[) can be replaced by 
|H^ A , where H^a = d^h u \ + d\h^ + d^hx^ is the field-strength tensor of h^. The 
integration over h. pu in the resulting theory with the Lagrangian 



LL 



(3.6) 



is straightforward (see e.g. [23] for details) and yields 



Z^=exp' bMt) 



+ ^ y d 4 xdV£^(x)An(x-*^(z') |x 

y [dS^exp - (j-^-J J d 4 xd 4 x'S pu (x)D m (x - x')S flu (x')+ 



+ 



(^) 2 M fe y ^xdVS^An^-*')^^') 



(3.7) 
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where D m = mK\(m\x\) / (47r 2 |x| ) is the Yukawa propagator. The formal measure [dS^] 
here, which replaces the measure [cffi^] in eq. ( |3.2| ), implies a certain prescription of 
the summation over world sheets of closed strings. For the dilute-plasma model of closed 
strings, this measure will be found in Section 5. 

According to the first term on the r.h.s. of eq. Q3.7| ), the Yukawa part of the potential 
satisfies the Casimir-scaling law, since M 2 k = C k =~ I ^ v fc) g3|. Further, in the leading 



semiclassical approximation we are considering in this and next sections, the integral over 
small-sized closed strings can be disregarded, and the remaining x E^-interaction 
produces the Casimir scaling also for the confining part of the potential. In fact, extracting 



from this interaction the string tension according to the respective general formula [33] 
(cf. also |25|]), we obtain a k = C^ct, where a = 4TrNi] 2 lnK. Note that, for o\ to be 
A^-independent as the quark-antiquark string tension in QCD, we should have 

V ~ -J^f= = 0{N- l > 2 ). (3.8) 
V-/V In k 

In what follows, we will address corrections to the Casimir scaling, that appear from the 
deviation from the London limit, as well as corrections produced by closed dual strings. 

4. Corrections due to the deviation from the London limit 



Let us consider the Lagrangian (|2.4j ) with £* = 0. Introducing the h^-field in the same 
way as in the London limit, we obtain the following expression: 

where the tensor = ^0- ^ fiQj'Qi is apparently symmetric. Performing the Gaussian 

i 

integration over B^, we arrive at the following substitution: 

iB^V + + K — > ^ d "Kv) " ^) W^JO- ( 41 ) 

It can be shown (see Appendix A for details) that det _1//2 [l + 3] , where 1 and 3 are 
the unit and the £ a ^-matrices, produces a renormalization of m and m#, which does not 
violate the London limit condition, ln«; > 1. Notice also that we have retained only the 
term linear in 3 on the r.h.s. of eq. ([4.1|). As we will see below, this linear term eventually 
produces a correction to a k , which we are looking for. Instead, the omitted 3 2 -term is 
shown in Appendix A to produce merely an inessential correction to m#, which is smaller 
than ran in the factor O ( — \= 



IN j 

Next, the <5 Q ^-term on the r.h.s. of eq. ([4.1|) is clearly the kinetic term of the Kalb- 



Ramond field, that is the first term on the r.h.s. of eq. (3.5). Instead, the £ a ^-term on the 
r.h.s. of eq. (|4.1| ) is the Higgs-inspired correction, which we will denote as C a/3 £ al3 , where 
C a P = — ■^i{d v h c ^ v ){d\h^ ^) . The Gaussian integration over ifi then leads to the following 
substitution: 
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(4.2) 

Analogously to eq. ( |2.3j ), the orthonormality of roots yields the following formula (cf. 
Ref. 0): 

E tfef ?7?i = 2 1} (<^ 75 + + ^ 5 ^ 7 ) • 

i 

Using this, we finally obtain the following Lagrangian: 

C k = -C aa + ±h%, - Im^VE^- 



7y 2 iV(iV + 1) 



d x D mH [x 



x') \c aa {x)C^{x') + 2C aP {x)C a ^{x') 



(4.3) 



To proceed, notice that, according to eq. (|3.3[), the physical meaning of the field is the 
monopole current, which couples to the dual gauge field B M and also possesses its own self- 
interaction terms. (Since for the massive vector field one always has c^B^ = 0, this current 
is automatically conserved.) The monopole current is defined in terms of the Kalb-Ramond 
field by means of eq. (|3.4[). Using this correspondence, it is straightforward to reformulate 



the action, corresponding to the obtained Lagrangian (4.3), in terms of monopole currents: 



j d 4 xC k = j cf^xk 2 + m 2 J d A xd A x'~k^(x)Do(x — x')'k^(x')+ 



+igmV2M k J d 4 xd 4 x't fJiU (x)'k u (x')d : {lD (x - x') - ^ N (j^ + fj j d 4 xd 4 x' D mH (x - x')x 

x 'k^(x)k^(x)4(x')4(x') + 2k«(x)kP(x)K{x')kP{x')] , (4.4) 

where the terms are presented in the same order as they stand in eq. (|4.3|) . The first 

and second terms on the r.h.s. of this equation are clearly the mass term of the current 

and the Coulomb interaction between currents, respectively. The third term has a form 

similar to the Gauss linking number of a closed surface S and a closed contour T, L(XJ, T) = 

J d xd 4 x''E fll ,(x)j l ,(x')d^DQ(x — x'),wheiej^(x) = § dx fl (T)8(x — x(t)). Our formula differs 

r 

from the Gauss one in two respects: first, the surface is open, and second, instead of the 
classical current of a point-like particle, localized along a closed trajectory, we have the 
quantum field, gmy^M^k^, distributed over the whole space-time. 

Let us now estimate the Higgs-inspired correction, given by the last term on the r.h.s. 
of eq. ( [4.4[) . Notice that, according to eq. ([3.8]) , m depends on N as m = g m r]\fN ~ 
V^V • • VN = v 7 ^, and mu = ton, where k has been naturally chosen iV-independent 
(cf. the first paragraph of Section 3). Therefore, mn grows with N as 
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m H = 0(N 1 ' 2 ). (4.5) 

The heaviness of the Higgs bosons, implied in both London and large- N limits, enables us 
to write, in the leading l/m//-approximation, the following expression for the correction 
under study: 



{r)m H ) 2 N(N + 1] 



J d*x [(g) 2 + 2(k,k,) 2 ] . (4.6) 



After the variation of the total action with respect to fc^, the color structure of the saddle- 
point equation thus obtained prescribes to seek k^ in the form M/c^. The equation for k^ 
then reads 

i 1 ~ (rjrn H yN(N + l) k t) *" + m * / ~ = 
= -M| J d *x% u (x')dZD (x - x'). 

The leading-order part of this equation, without the /c 2 -term, can easily be converted 
into the differential form by acting with d 2 onto both its sides. The resulting equation 
(d 2 — m 2 )k fl = % -2j=-d L ,Y>p U leads to the following leading-order saddle-point expression for 

kfj,: k fl (x) = J d 4 x'D m (x — x')d u T,^ u (x'), which should further be substituted into 

eq. Q, 



d 4 x(kl) 2 . (4.7) 



{r)m H ) 2 N{N + 1) 

Among all the terms contained here in A; 2 , only the surface x surface one, 

_(gmf I ^ Xl ^ X2 ^ Xl) ^ X2) ^ Dm ^ _ Xl ) 9 ^D m (x - x 2 ), 

produces the desired correction to the string tension. Indeed, the integral structure of the 
respective part of the correction (^.7|), 

3(Cfc) 2 (gm) 4 f 

da^{xi)da^(x 2 )daxp{xz)daxp{xi) • J, (4.8) 



8N(N + l)(r)m H y 

4 

where J = d* 1 d% 2 <92 3 d^ 4 f d 4 x Y\ D m (x — xi), is the same as appears in the 3-d SU(iV) 

i=i 



Georgi-Glashow model due to the non-diluteness of the monopole plasma [13]. In that 
paper, it was shown that this structure does produce a correction to the string tension. 
It can be shown that, similarly to the 3-d SU(iV) Georgi-Glashow model, this correction 
behaves with N as (C k ) 2 /N. 

With the account for the obtained correction, the iV-ality dependence of the ratio of 
string tensions is given by the formula 
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a k _ k(N - k) 
Vi ~ N - 1 



N 

where the coefficient a(N) ~ A -1 and therefore, at fixed k, the whole correction vanishes 
in the large-iV limit. As well as the Casimir-scaling term, the obtained leading Higgs- 
inspired correction is apparently invariant under the interchange of quarks and antiquarks, 
k «-> (N-k). 

5. Corrections due to closed strings 

To study the grand canonical ensemble of closed strings, it is necessary to replace 

in eq. (|3.6| ) (with set for a while equal to zero) by Yl n a$%, where the n a stand 

a 

for winding numbers. It is known [31, |2^] that one may restrict oneself to closed strings 



possessing the minimal winding numbers, n a = ±1. That is merely because the energy of 
a single closed string is a quadratic function of its flux, owing to which it is energetically 
favorable for the vacuum to maintain two closed strings of a unit flux, rather than one 
string of the double flux. 

Then, taking into account that the plasma of closed strings is dilute, one can perform 
the summation over the grand canonical ensemble of these objects, which modifies the 
Lagrangian fl3.6| ), with T,^ = 0, as follows: 

Vca*. = i^H^ + ig„ - 2Ccos (f^) • (5.1) 

Here A = \J Lja? is a UV momentum cut-off with L and a denoting the characteristic 
distances between closed strings and their typical sizes, respectively. In the dilute-plasma 
approximation under study, a -C L and A 3> a -1 . Next, Q oc e~ s ° stands for the fugacity 
(Boltzmann factor) of a single string, which has the dimension (mass) 4 , with So denoting 
the action of a single string, So ~ (J\a 2 . Finally, it is assumed that closed strings are not 
too small, namely a > O I J , so that So 3> 1, and the mean density of the plasma, 2£, 

\ 9 m J 

is exponentially small, i.e. the plasma is dilute. 

Note also that, because of the Debye screening of the dual vector boson in the plasma 



of closed strings, its mass increases. This is clearly seen from eq. (5J), by the increase of 
the mass of the Kalb-Ramond field, which represents this boson: 

m 2 — M 2 = m 2 (l + . (5.2) 



4A 4 , 

To study corrections to the /c-string tension produced by closed strings, we will need to 
know correlation functions of these strings in the plasma. To obtain an expression for the 
generating functional of such correlation functions, one needs the theory to be formulated 
in terms of dynamical vorticity tensor currents. This can be done by recalling that, for 
closed strings: 



exp < — / d A x 



12m 2 ^ A ^4 ^ 
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VSnu exp 



d A xd 4l x'S fiU (x)D rn (x - x')S^(x') 



4 



j d A xh 



The Kalb-Ramond field can then be integrated out by solving the saddle-point equation 
stemming from the respective part of the Lagrangian, 



— 2£ cos 



gKA) ig h s 

4 A 2 / 4 n ^°^ - 



This yields the following expression for the partition function of the grand canonical en- 
semble of closed strings in terms of their vorticity tensor currents: 



gr. can. 



VSav exp 



/ ^ 4xrf VS Ml ,(x).D m (x - x )S fJiU (x') + y[S Mi/ ] 



where the potential ^[S^,,] reads 



(5.3) 



V[S 



d A x < A 2 1 | In 



A 2 , 
— S 



2C A /1+(^|S M ,| X2 



(Note that the operator, which describes the density of plasma at the point x, is A 2 |S At!/ ( 



(5.4) 



Comparing now eqs. (|5.3|) , ( |5.4j) with eq. fl3.7|) , we see that, in the dilute-plasma model of 
closed strings, the formal measure [cZS^] concretizes as VS^e^^^. 

Corrections to the string tension of the open world-sheet £ stem from the last term 
in the following expression for the partition function Q3.7|) : 



lnZ fc LL 



(gM k f 



j> dx^ j> dx'^D m (x - x 1 ) + ^- J d A xd A x'T,^(x)D m (x - x')Y^ u (x) 



c c 



In ( exp 



(^) 2 M fc J d A xd A x'S^(x)D m (x - x')^ u (x') 



where the average (...) over closed strings is now defined according to eqs. (|5.3j), (|5.4| ). By 
virtue of the cumulant expansion, this term can be written as 



71=1 



x / d A x 1 ---d A x n a^ Ul (x 1 )---a flnUn (x n )((S^ Ul (x 1 )--- S% Un (x n ))) 
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Here, a^ u (x) = J d A x'D m ( x — x')S^jy(x'), and ((■■■)) denotes a one-particle irreducible 
average (cumulant) of closed strings. Since the action corresponding to the partition func- 
tion ( |5.3|) contains only powers of S flLJ (x)S PjL/ (x'), cumulants of odd orders vanish, whereas 
a cumulant of an even order n has the form 

({s^M) ■ ■ ■ s Zu n M)) = 



■■■5° 



-in 



n flilVx,H2V2\ X l x 2) ' ' ' fn 



Xl p&2 



\ x n— 1 



Xr, 



+ permutations. 



Here, / w 

ui,fi 2 u 2 ( x i ~ x 2) — £ fj,ivi\ip £ ^2^2^2P^xl^\l ^{ x i ~ x 2) with T> standing for a function 
whose concrete form for the case of a very dilute plasma will be made clear in a moment. 
As for the Lorentz structure of the /-tensors, it stems from the condition of closeness of 
strings, d^S^u = 0. The color structure of the cumulant produces, for n = 21, the factor 
(M|)' = (Ck) 1 ■ The iV-ality dependence of is therefore defined by the following formula: 



(5.5) 



1=1 



where cr^ 1 ' are (iV-dependent) coefficients of dimension [mass] 2 . Note that the term with 
I = 2 in this equation produces a correction to the string tension, which has the same 
(Cfc) 2 -dependence as the leading Higgs-inspired correction found in the previous section. 

As an example, let us finally present the lowest non-trivial two-point correlation func- 
tion of closed strings, which can be derived in the approximation when the plasma is very 
dilute, i.e. its density is even lower than the (already exponentially small) mean one, 2£. 
In that case, the potential ( p\4| ) becomes a quadratic functional. Including the source term, 
J d iJpSp, into the square brackets on the r.h.s. of eq. (^), we obtain for this Gaussian 
integral: 



-2^gr. can. [J/zi/] 



^gr. can. [0] 



VSav exp 



(X)V dAxdix 'SA X ) D rni x - x ')SA X ') + 



+ I d*x( -2C + ^S 2 , + J M „S Mi , 



exp 



J d A xd i x'3 tMU (x)Q{x - y)J^(x') 



where Q{x) = jr(d — m )Dm(x), and the mass M is defined by eq. ( |5.2| ). Imposing the 



condition d^S^ = 0, one can further, similarly to Ref. [29|, derive the desired two-point 
correlation function (string propagator): 



C 



jd^(d 2 -m 2 )[D (x)-D M (x)}. 



(5.6) 



(MA 2 ) 2 

It is known that, in this Gaussian approximation, all higher cumulants vanish, i.e. the terms 
with I > 2 in eq. ( |5.5[) are absent. Therefore, in this very dilute plasma approximation, the 
Casimir-scaling law is preserved, since the terms that violate it vanish. 
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6. Summary 



In this paper, we have explored the fc-string tension spectrum in the SU(iV)-inspired 4-d 
dual Abelian-Higgs-type theory. We have first considered the London limit of this theory 
and demonstrated that, in the leading semi-classical approximation when the small-sized 
closed dual strings are completely disregarded, the /c-string tension obeys the Casimir- 
scaling law. In the same approximation, when closed strings are disregarded, we have 
further explored the leading correction to the Casimir scaling emerging due to the deviation 
from the London limit, i.e. due to the finiteness of the masses of the Higgs bosons. This 
correction turns out to have the same iV-alrty dependence as the correction, which one 
finds in the 3-d SU(iV) Georgi-Glashow model when the non-diluteness of the monopole 
plasma is taken into account. We have then addressed another type of corrections to the 
Casimir scaling, which emerge in the London limit when one accounts for the dilute plasma 
of closed dual strings. In the leading low-density approximation, i.e. when the plasma is 
very dilute, the respective correction is shown not to violate the Casimir scaling. Instead, 
the correction of the next order in the non-diluteness has the same iV-ality dependence 
as the above-mentioned correction emerging without closed strings in the vicinity of the 
London limit. Finally, we have analyzed the corrections that appear in higher orders in the 
non-diluteness of the plasma of closed strings. Interestingly, the 1/N dependence of the 
ratio (Jk/cri does not satisfy the counting rules that were spelled out for S\J(N) Yang-Mills 



theories in Ref. [12]. However, one has to be careful in comparing our results to the full 



non-Abelian gauge theory. In our model, off-diagonal degrees of freedom are disregarded, 
the only remnant of the non-Abelian structure being the quantization condition for the 
magnetic charge. The discrepancy between the results presented here (and in Ref. []i~3| l) 
and those obtained in Ref. [^] suggests that a purely Abelian description of the SU(iV) 
vacuum is not adequate to catch the full dynamics of the non-Abelian gauge theory, in 
agreement with the conclusions in Refs. [05| , 36]. Further work along these lines is needed 
to clarify this issue. 
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A. A few technical points 



Let us first evaluate the effect produced by the determinant when integrating over the 
B^-field in eq. Q4.1|) . One has 



det" 1/2 [1 + H] ~ exp 



•-TVS + -Tr~ 2 
2 4 



(A.l) 
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where "Tr" includes both the trace over the indices a, (3 and the trace in the coordinate 
space. Next terms have been omitted since they describe self-interactions of the y^-field, 
rather than the renormalization of m and run , which will be shown to be produced by the 
two retained terms. The first term on the r.h.s. of eq. ( A.l| ) modifies eq. ( |4.2| ) as follows: 



m 



(ja/3 _j_ zm H 



d x Djyijj (^x 



Cr"(x') + ^JLp s 



\x=0 



zmjr. Note 



where the regularization parameter z is defined by the relation 6^(3 
that, due to eq. Q4.5Q and the fact that S^(x)\ x=Q is iV- independent, 2 should scale with 
iV as 0{N~ 2 ). Apart from this requirement, z can be chosen at will. Up to an inessential 
constant addendum, the correction to the action thus reads: 



2zml 



§ E / d^x'D mH {x - x 1 ) \c^(x)f< s + CT 5 (x')6^ 



2zm 2 H 
' Nr] 2 



d A xC c 



Comparing this result with the first term on the r.h.s. of eq. ( [4.3[ ), we arrive at the following 
renormalization of m~ 2 : 



1 



1 

m 



2zm 2 H 
> \ ^ Nr] 2 



(A.2) 



2zrn 2 H 



0(N 1 ). Let us now consider the 



where the correction vanishes at large N as N ^ 
second term on the r.h.s. of eq. ( |A.1| ). It modifies the mass term on the l.h.s. of eq. (| 

as 



m 



H ,„2 
Pi 



Azm H „<x3„<x3 



(N V ) 



2 1i '<i 1j 1j 



The tensor qfq^qfqj should be proportional to 5ij, which is the only tensor symmetric in 
indices and the proportionality coefficient is 1. The squared Higgs mass therefore 

renormalizes as 



m 2 H 



m 2 H 



Azm 2 H 



(Nr}f 



(A.3) 



The obtained correction vanishes at large N: 



4zm 2 r, 



0(N~ 2 ), so that < 1 in the 



WW ~ v ; ' WW 
large- N limit. A more accurate condition on z, which should hold at any N, can be imposed 

by using the inequality of the London limit, eq. ( |3.1[ ). It yields: 
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4zm 2 H _ 16Xz 16z (2vreA0 2 _ (8vre) 2 z 
> Jmif ~ ~W > W A ~ A : 

i.e. at a given A, z should be chosen such that z <C ^ 8 \ 2 • The squared Landau-Ginzburg 



parameter then renormalizes according to eqs. ( |A.2j ), ( |A.3| ) as 



2 2 
K — y K 



1 + i 1 - 1)] = + ° (Ar_1) • t 1 + 0(iV " 1 )] > • 



Therefore, upon renormalization, k remains large, and the respective correction vanishes 
in the large- N limit. 



Let us now evaluate the omitted 3 -term on the r.h.s. of eq. (4.1). This term reads 



k^^kl = —k«k; m q?q?q?q], (A.4) 

with the monopole current k M defined by eq. ([O]). Since the tensors k^k"^, tpifj are 
symmetric in indices (a, 7) and respectively, it is natural to impose the following 

Ansatz: q^q^q^qj = M5 ai 5ij. The proportionality coefficient J\f can readily be found: J\f = 

8k?, 



and the term (A..4) therefore reads (^jv) a (Ar-i) ^' Furthermore, the characteristic 



N- 

amplitude of the k^-field can be estimated by noticing that the configuration of this field, 
which dominates in the partition function, is the one at which each of the first two terms 
on the r.h.s. of eq. ( ^4.4| ) is of the order of unity. When applied to the mass term, this 
requirement yields L ~ , where L and |k^| are the characteristic wavelength and 

the amplitude of the k„-field, respectively. Applying further the same requirement to the 



Coulomb interaction of monopole currents, the second term on the r.h.s. of eq. (4.4), we 
have m^k^l 2 X 6 ~ 1. Substituting here the above estimate for L 2 , we obtain the desired 
estimate for the characteristic value of the amplitude: |k^| ~ m 2 . This leads to the 
following estimate for the magnitude of the term ( |A.4[ ): fajvW]V-l) ^i' This term therefore 
produces a small positive correction to run, whose magnitude with respect to m# can be 
estimated as: 

m 2 1 f 1 

O 



rjNy/N -lmfl \kVn , 
B. More on Casimir scaling in the 3-d SU(iV) Georgi-Glashow model 



It has been proved in [13j that, in the 3-d SU(iV) Georgi-Glashow model, Casimir scaling 
holds for an arbitrarily shaped surface (i.e. the world sheet of a /c-string), provided that 
the density of monopole plasma is much lower than the mean one. The purpose of this 
Appendix is to show that, in case of a flat surface, Casimir scaling in this model is an 
exact result, not requiring the condition that the density of the monopole plasma is much 
lower than the mean one. To this end, let us consider the confining part of the k-th power 
of the fundamental Wilson loop, i.e. the part produced by monopoles. It has the form 

TV 

a^(C)) mon = E W ai ,..., ak (C), where 

( a lt ...,a h =l \ 

Vwith possible coincidences/ 
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(C) 



Zn 



VB^ie^xduBx) / Plexp^ / d 6 x 



9m 



2 K+ 



+ 47riM 



(n) 



dcTfjB^ > , 



(B.l) 



S(C) 



where ^ mon is the same functional integral, but with the last term [which describes the flux 
of the magnetic field through an arbitrary surface S(C)] set equal to zero. The constraint 
tfiuxduBx = imposes the fact that free photons, inessential for confinement, are not taken 
into account. The dimensionalities of the magnetic coupling constant, g m , dual photon field, 
1, and the (exponentially small) monopole fugacity £ are [mass] -1 / 2 , [mass] 1 / 2 , and [mass] 3 , 
respectively; for more details on the model and eq. ( |B.1| ) at k = 1 see [13, |3"3} ]. In eq. ( |B.1| ), 



M 



^ again denotes the sum ^ m a t , where some n indices out of k can now coincide. 

i=l 



For the contour C located in the (x,y)-plane, the saddle-point equations stemming 
from eq. ( [B.lD read 



ig m l' + g 2 m B-4mM k n) 5(z) = 0, 



(B.2) 



iB' -2(Y,<li MdmCiil) = 0, (B.3) 

i 

where ' = d/dz, and the natural Ansatz B M = 5^B(z), 1 = \{z) has been adopted. Next, 
as it follows from eq. flB^) , B cx iM.[ n) , therefore W ai ,..., ak {C) -► e"^"') ct|S(C)I at 



|S(C)| — > oo, where the string tension a is ^-independent. Since ^M^J = Cfc + w '^"" , in 
the limit of asymptotically large areas |S(C)| of interest, we arrive at a Feynman-Kac-type 
formula, 



,2_ 



k 



(W k (C)) mon = e- c -| 2 (c)| ^ Cne -^^(C)| (B . 4) 



n=l 



with some positive coefficients c n ; therefore, only the case n = is relevant in eq. (B.l) 
(cf. ref. [|DJ). We should, thus, solve the system of eqs. ( |B.2| ), ( |B.3| ) with M^, which 
coincide with from the main text. Setting B(z) = ~M.i-B(z), l(z) = MfcZ(z), we see that 
eq. ( p.2| ) takes the same form as in the fundamental case, namely 

ig m l' + g m B = 4m5(z), (B.5) 

whereas to handle eq. ( |B.3| ), one should notice that any root vector is a difference of two 
weight vectors, = q a b = m a — m^. In particular, positive roots, we are dealing with, are 
those with b < a, therefore eq. ( |B.3[) takes the form 



^2 M fc (m Q - m 6 ) sin [g m M k (m a - m b )l] = 7^C k B' . 

b<a ^ 
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Using the symmetry of the expression under the sum on the l.h.s. with respect to a <-> b, 
we can rewrite this equation as 



N 

M fe m a [sm(g m M k m a l) cos(g m M k m b l) - cos(g m M k m a l) sin(g m M fc m 6 Z)] = —C k B'. 

a,b=l * 

(B.6) 

We should further perform the four sums on the l.h.s.; let us begin with the first one, 

N 

^2 (Mfcm a ) sin(<7 m Mfcm a Z). Appartently, there are k terms in this sum, for which m a 

a=l 

coincides with some of the weight vectors, which enter M^. Using the relation m a mft = 
(5 a b — N~ 1 )/2, we have for such terms M. k m a = — (k — 1)^ = t^\T- For the other 
(N — k) terms in the sum, m a does not coincide with any weight vector in M^, hence 
Mfcm a = — 2^ for such terms. We, therefore, obtain 



2^(M fc m a ) sm(j m M t m a i) = k ■ ^ N sm [9m ^ I I 

a=l ^ ' 



2iV < +^-^-2iV Sin l^ 



fc A . f N-k. 
sm l g m -^-:l + sin g m . — t 



2N J V 2Ar 

Remarkably, already this expression alone is manifestly invariant under <-> (iV — fc). In 
the same way, we obtain for the three other sums the following expressions: 



^ cos(# m M fc m fc /) = k cos I # m 
6=1 ^ 

N 

yj(M fc m a ) cos(ff m M fc m a Q = C fe 



N-k 
2N 



I ) +(JV- A;) cos I 5m^Z ) , 



a=l 



cos 5. 



N-k 
' 2N 



I - cos g. 



'2JV 



"2iV 



iV 



M fc m b /) = sin I g, 



N-k 

' 2N 



I) -(N-k)sm[g m —l) . 



k 



The l.h.s. of eq. (|B.6|) then takes the form C^iV sin and the whole equation becomes 



B' + 2i(Nsm^ = 0, 



(B.7) 



that also coincides with that of the fundamental case p4| ]. The solution to the system of 
equations ( B.5[ ), ( p.7| ) reads 



„, . .8m D e - m °W „ . 8 

B(z) = 1 — k — ; — r , l(z) = — sen z ■ arctan e 

g 2 m l + e- 2m ^l^l' g m 1 



-m D \z\ 



(B. 



where m_o = g m V^C 1S the Debye mass of the dual photon 1. [This mass is visible in 
eq. (|B.1|) with C = 0, where the B M -field is integrated out to produce the kinetic term of 
the 1-field, and cosine is expanded up to the quadratic term.] Inserting the so-obtained 
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field B(z) into eq. (|B,1| ), we obtain a = cgy/N£, where g = A^jg m . The proportionality 
coefficient c here depends on the range of average of B(z), since, according to eqs. flB,g| ), 
the string is exponentially thick, \z\ < wi^. F° r instance, if we simply choose the value 
-B(O), then c = 4. [For comparison, the value of c one obtains in the case when the density 
of the monopole plasma is much lower than the mean one, £N(N — 1), is |p3 tt, that can 



be shown to approximately correspond to the following range of average: \z\ < 
The obtained string tension is manifestly ^-independent. 

Equation ( |B.4j) then leads to the conclusion that, for a flat surface, Casimir scaling 
holds in the full sine-Gordon theory describing the monopole plasma in the 3-d Georgi- 
Glashow model. In another words, for a flat surface, Casimir scaling holds also at monopole 
densities close to the mean one, £iV(iV — 1), rather than only at the densities much smaller 
than (N(N — 1), as it takes place for a non-flat surface (l3j . It is finally reasonable to 
have some feeling on how much a surface should deviate from a flat one in order that the 
Casimir scaling starts violating, if the density of monopoles is the mean one, £N(N — 1). 
To this end, let us consider a straight string (apparently corresponding to a flat surface) 
of a minimal possible length, m^ 1 . According to the first of eqs. (|B.8| ), the field B(z) at 
the end points of such a string decreases above and below the surface also at the distance 
nip 1 . Therefore, if we start bending the string such that it forms a piece of a circle, the 
two solutions overlap with each other if the radius of this circle is < m7) ■ The critical 
situation when this radius is equal to m^ 1 corresponds to the distance between the end 
points of the string equal to 2m^ 1 sin | . This is the minimal distance which should hold 
between two points, separated by the distance m^ 1 along the string, at which the surface 
can still be considered as flat from the point of view of the Casimir scaling. 
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